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Abstract — This  paper  presents  an  interior-point  method  (IPM) 
based  quadratic  programming  (QP)  solver  for  the  solution  of 
optimal  control  problem  in  linear  model  predictive  control 
(MPC).  LU  factorization  is  used  to  solve  the  system  of  linear 
equations  efficiently  at  each  iteration  of  IPM,  which  renders 
faster  execution  of  QP  solver.  The  controller  requires  internal 
states  of  the  system.  To  address  this  issue,  a  Proportional 
Integral  Observer  (PIO)  is  designed,  which  estimates  the  state 
vector,  as  well  as  the  uncertainties  in  an  integrated  manner. 
MPC  uses  the  states  estimated  by  PIO,  and  the  effect  of 
uncertainty  is  compensated  by  augmenting  MPC  with  PIO- 
estimated  uncertainties  and  external  disturbances.  The 
approach  is  demonstrated  practically  by  applying  MPC  to  QET 
DC  servomotor  for  position  control  application.  The  proposed 
method  is  compared  with  classical  control  strategy-PID 
control. 

Index  Terms — Model  Predictive  Control,  Interior-Point 
Method,  LU  factorization,  DC  Servomotor,  Proportional 
Integral  Observer. 

I.  Introduction 

This  paper  considers  Interior-point  algorithm  used  on- 
line with  Linear  Model  Predictive  Control  (MPC).  Linear  MPC 
assumes  a  linear  system  model,  linear  inequality  constraints 
and  a  convex  quadratic  cost  function  [1].  MPC  can  be 
formulated  as  a  quadratic  programming  (QP)  problem  and 
solved  at  each  sampling  interval.  It  thus  has  the  natural  ability 
to  handle  physical  constraints  arising  in  industrial 
applications  [2]. 

At  each  sampling  instance  MPC  solves  an  online  QP 
optimization  problem,  computes  the  sequence  of  optimal 
current  and  future  control  inputs  by  minimizing  the  difference 
between  set-points  and  future  outputs  predicted  from  a  given 
plant  model  over  a  finite  horizon  in  forward  time.  Then,  only 
the  current  optimal  input  is  applied  to  the  plant.  The  updated 
plant  information  is  used  to  formulate  and  solve  a  new  optimal 
control  problem  at  the  next  sampling  instance.  This  procedure 
is  repeated  at  the  each  sampling  instance  and  the  concept  is 
called  as  receding  horizon  control  MPC  strategy.  Since  this 
quadratic  program  can  be  large  depending  upon  control 
problem,  MPC  requires  long  computation  times  at  each 
sampling  instants,  therefore  it  is  usually  restricted  to  systems 
with  slow  dynamics  and  large  sampling  intervals,  such  as 
chemical  processes  [3] .  The  ability  to  solve  the  QP  problem 
online  become  critical  when  applying  MPC  to  systems  with 
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fast  response  time  and/or  embedded  applications  where 
computational  resource  may  be  limited  [4] . 

Recently,  many  reports  in  the  literature  address  applying 
MPC  to  control  applications  with  short  sampling  intervals, 
by  adapting  fast  optimization  algorithms.  These  algorithms 
solve  the  resulting  QP  by  exploiting  the  special  structure  of 
the  control  problem  the  MPC  sub  problem.  IPM  approaches 
the  solution  of  Karush-Kuhn -Tucker  (KKT)  equations  by 
successive  decent  steps.  Each  decent  step  is  Newton's  like 
step  and  the  solution  is  obtained  by  solving  system  of  linear 
equations  using  appropriate  numerical  methods  in  order  to 
determine  search  direction.  Matrix  factorization  approach 
provides  a  means  to  simplify  the  computation  involved  in 
linear  solver.  Different  matrix  factorization  methods  are  used 
such  as  Gauss  elimination,  QR,  LU,  and  Cholesky  [3],  [5],  [6]. 

Model  predictive  control  has  received  considerable 
attention  driven  largely  by  its  ability  to  handle  hard 
constraints.  An  inherent  problem  is  that  model  predictive 
control  normally  requires  full  knowledge  of  the  internal  states. 
In  the  physical  implementation  of  control  strategy  based  on 
the  system  states,  if  the  states  are  not  available  for 
measurement  then  state  feedback  will  not  be  the 
implementable. 

The  problem  of  robust  control  system  design  with  an 
observer  based  on  the  knowledge  of  input  and  output  of  the 
plant  has  attracted  many  researchers  now-a-days  [7],  [8],  [9], 
[10],  [11],  [12],  [20] .  The  observer  proves  to  be  useful  in  not 
only  system  estimation  and  regulation  but  also  identifying 
failures  in  dynamic  systems.  This  requirement  of  obtaining 
the  estimates  of  uncertainty,  as  well  as  state  vector  in  an 
integrated  manner  is  fulfilled  by  PIO.  PIO  approach  basically 
works  when  the  system  uncertainty  varies  slowly  with  respect 
to  time.  It  is  an  observer  in  which  an  additional  term,  which  is 
proportional  to  the  integral  of  the  output  estimation  error,  is 
added  in  order  to  estimate  the  uncertainties  and  improve  the 
system  performance  [7],  [8],  [9]. 

In  this  paper,  we  present  an  interior-point  algorithm  to 
solve  MPC  problem,  which  utilizes  Mehrotra's  predictor- 
corrector  algorithm  [13],  and  the  linear  system  at  each  IPM 
iteration  is  solved  by  LU  factorization  based  Linear  Solver. 
PIO  is  designed  for  state  and  uncertainties  estimation.  The 
remainder  of  the  paper  is  organized  as  follows:  section  II 
describes  MPC  problem  formulation.  Section  III  presents  the 
interior-point  method.  In  section  IV,  the  PIO  theory  is  briefly 
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reviewed  and  applied  for  state  and  uncertainty  estimation. 
Mathematical  modeling  of  DC  servomotor  is  derived  in 
Section  V.  IPM  QP  solver  results,  the  effectiveness  of  observer 
based  MPC  scheme  with  simulation  results  is  presented  in 
section  VI.  MPC  implementation  on  Quanser  QET  DC 
servomotor  is  presented  in  section  VII. 

II.  Mpc  Problem  Formulation 

A  discrete  time  linear  time-invariant  model  of  a  system  in 
a  state  space  form  is  given  as, 


(1) 


x(t  +  1)  =  Ax(t)  +  Bu(t)] 
y(t)  =  Cx(t)         J' 

where,    y(t)  e  K.     are  output  vector,    «(f)  e  ^ 
are  input  vector  and    x(t)  e  ]f^n  are  internal  states  vector. 

A  and  B  are  system  parameter  matrices.  Given  a  predicted 
input  sequence,  the  corresponding  sequence  of  state 
predictions  is  generated  by  simulating  the  model  forward 

over  the  prediction  horizon  (Np )  intervals. 

x(t  +  2)  =  Ax(t  +  l)  +  Bu(t  +  l)        ) 
x(t  +  2)  =  A2x(t)  +  ABuit)  +  Bu(t  + 1) 


x(t  +  NP)  =  AN'x(t)  +  £,:;  (AN--l-jB)u(j). 
The  prediction  model  is  given  by, 

y(t  +  NP)  =  Cx(t  +  NP). 
Putting  the  value  of  x(t  +  NP )  in  (4),  we  get, 

y(t  +  NP)  =  C[AN"X(0)  +  ^(AN^-JB)u(j)]. 

In  matrix  form  (5)  is  written  as, 

Y  =  ®x(0)  +  FU. 
with, 

Y  =  [y(t  +  l)     y(t  +  2)    ...     y(t  +  NP)J, 
U  =  [u(t)     u(t+T)    ...     y(t  +  Nu-l)J, 


(2) 
(3) 

(4) 


(5) 
(6) 


O 


CA 
CA2 


CA1 


CB 

0 

0 

r  = 

CAB 

CB 

0 
0 

CAN"lB 

CAN-- 

lB     ■ 

•     CAN'-N»B 
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where,  N u  is  control  horizon.  The  goal  of  the  controller  is  to 
make  the  difference  between  the  output,  y(t  +  Np  )  and  the 

reference,  yrei  (t  +  Np  )  as  small  as  possible.  This  can  be 

defined  by  using  a  standard  least  squares  problem. 
The  objective  function  is  defined  as, 


J=^  2N + ^  -  y<^ + ')|L + ~  Elh^ + C-    (7) 

where,  Au(t  +  i)=u(t  +  i)  —  u(t  +  i  —  1) 

subjected  to  linear  inequality  constraints  on  system  inputs, 

u       <  u(t)  <  u 

where,    y ' re*  is  the  set-point  and  while  Q  and  ^  are  the 

positive  definite  weight  matrices  [14].  To  get  optimal  solution, 
the  MPC  problem  can  be  formulated  as  a  standard  quadratic 
programming  (QP)  problem  as, 

/([/)  =  \~UTHU  +  /rC/|subjected  to  AU  -b  <  0. (8) 

where,  //  is  (IN   x  IN  )  Hessian  matrix,  f  is 
(INU  x  1)  column  vector  [14],  [16]. 

III.  Interior-Point  Method 
Consider  the  following  standard  QP  problem, 

mj11  ±xTQx+CTx subjected  to  Ax  - b  =  0, Lx  -  k  <  0. 

(9) 
Except  equality  constraints,  (9)  resembles  (8).  The  KKT 
Conditions  for  optimality  are  given  by, 


Qx  +  ATy  +  LTz  =  -C 
Ax-b=0 


Lx  +  s-k  =0 
0 


T 

z  s 


z,s>0 


(10) 


where,  y  and  z  represents  the  Lagrange  multipliers  for 
equality  and  inequality  constraints  respectively.  The  slack 
variable  5  is  introduced  for  converting  inequality  into 
equality.  By  applying  Newton's  method,  search  direction 

(Ax,Ay,Az,  As)  is  obtained  as  shown  in  (11),  where 
J(x, y, z, s)  denotes  Jacobian  of  f(x,y,Z,s)  , 


J(x,y,z,s) 
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Ax 

1 

Ay 

h 

Az 

r3 

As 

_r4_ 

(11) 
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Considering  only  inequality  constraints,  augmented  form  of 
original  system  is  given  as, 


(12) 


Q 

LT 

Ax 

ri 

L 

-ZlS 

_Az_ 

_r3-Z~lr4_ 

As  =  Z-l(r4-SAz). 


(13) 


We  use  Mehrotra's  Predictor-corrector  algorithm  [13]  to  solve 

(12)  and  (13)  for  (Ax,  Az,  As).  Predictor  or  affine  and 

corrector-centering  steps  uses  (r\ ,  r3 ,  r4  )  values  from  (14) 
and  (15)  respectively  to  form  the  right-hand  side  of  ( 1 2)  and 

(13)  which    in    turn    finds  (Ax ajf '  ,Az"ff '  ,Asajf  )  and 
(Axcc ,  Azcc,Ascc )  [1],  [15],  [17],  [18],  [19]. 


(14) 


Qx  +  C  + 

L  z 

Lx  +  s  - 

-k 

ZSe 

0 
0 


(15) 


diag  (Azaff  )diag  (Asaff  )e  -  o/ue 
where,  <J  is  centering  parameter  and  JU  is  complementarity 
measure  or  duality  gap.  Thus,  solution  of  linear  system  in 
(13)  is  to  be  calculated  twice  at  each  iteration  of  IPM,  which 
puts  major  computational  load  on  overall  performance  of  IPM, 
and  eventually  on  MPC.  Linear  system  in  (13)  uses  two 
different  right-hand  side  components  but  the  same  coefficient 
matrix,  so  just  one  factorization  of  this  matrix  is  required  per 
iteration.  Such  system  can  be  solved  using  Gauss  elimination, 
QR,  LU,  Cholesky  factorization  methods. 

The  coefficient  matrix  in  a  linear  system  (13)  involved  in 
IPM  is  symmetric  and  indefinite  so  conventional  and  more 
numerically  stable  Cholesky  factorization  cannot  be  used  to 
solve  the  linear  system.  Compared  to  computational  cost  of 
QR  decomposition  (2/3n3)  and  the  problem  of  associated 
pivoting  with  ill  condition  of  Gauss  elimination,  LU 
factorization  proves  to  be  more  effective  and  accurate  to  solve 
system  of  linear  equations  [5],  [6].  We  propose  LU 
factorization  based  Linear  Solver  for  the  same.  The  pseudo 
code  of  LU  factorization  Linear  Solver  for  (13)  is  given  as 
follows: 

Pseudo  Code-LU  factorization 
A  =  |An     A12;     A21     A22J; 

where,  Au  =  Q;An  =  LT;A22  =  L;A22 

B  =  [Bn;     B22], 

where,  Bn  =  rx,B22  =  (r3-Z^r4); 

©2013  ACEEE 
DOI:01.UCSI.4.1.1064 


-z  ls 


A  =  LU; 

Solve  LY  =  B;  //forward  substitution 
UX  =  Y;  //backward  substitution 
X  =[Xl;     X2\; 


Interior-point  method(IPM)  algorithm 

Initialization  and  input 

X0  =  [0    0     ...     Of; 

Z0  =  [l    1    ...    lF;S0  =  [l    1    ...    lf;ZQ,SQyO; 

e  =  [l    1    ...     lf;a=0.25 

Input  Q,C,L,k  matrices  from  given  QP  problem 

o     o 
/U  = ;  m-  no.  of  constraints 

m 
Start  the  loop  and  terminate  if  stopping  criterion  are  satisfied 
Predictor  or  affine  step 

S  =  diag(S0);  Z  =  diag(Z0); 

ruff  =  QXo  +C  +  L  Z0;     r3aff  =LXQ  +  S{)  - k; 

r4ajf  =  ZSe; 

Solve    (13),    (12)    by    linear    solver    and    compute 

Axaff  ,Azaff  ,Asaff 

Compute  centring  parameter  <T 


Mali 


a 


(Z0+aAzaff)T(S0+aAsaff) 


in 


f  ..     A 


Vaj 


v  M  j 


Corrector  and  centring  step  to  obtain  search  direction 

Zaff=diag(Azaff);  Saff  =diag(Asaff); 

ricc  =  0;        r3cc  =  °;        rAaff   =  Zaff  S off  e  ~  OfW, 

Solve    (13),    (12)    by    linear    solver    and    compute 

Axcc,Azcc,Ascc 


Update  (Ax,Az,As)  as 

(Ax,  Az,  Aj)  =  ( Ax"# ,  Az"ff ,  As"ff )  +  (Axcc  ,Azcc ,  Ascc ) 

Update  XQ,Z0,S0as 

XQ  =  XQ+aAx;    ZQ=ZQ+aAz;    SQ=SQ+aAs; 
End  of  loop 
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IV.  Observer  Design 

The  implementation  of  MPC  requires  the  measurements 
of  all  the  states.  As  only  an  encoder  is  used  for  measurement 
of  the  motor  shaft  position  in  practice.  PIO  is  designed  in  this 
section  for  present  application. 

Consider  a  linear  time-invariant  system 


x(t  + 1)  =  Ax{t)  +  Bu(t)  +  L(y(t)  -  y(t))  +  Beit) 
y(t)  =  Cx(t) 

e(t  +  \)  =  e(t)  +  M(y(t)-y(t)) 


■  (16) 


where  x(t)  is  the  estimated  state  vector  and  e  (t )  is  the 

estimation  of  uncertainties  of  p-dimensional  vector.  L 
proportional  and  M  integral  observer  gain  matrices  of 
appropriate  dimension. 

The  system  described  by  Eq.  ( 1 6)  is  said  to  be  a  full-order 
PI  observer  for  system  (1)  if  and  only  if 


limx(0  =  0l 

»-»00  \ 

lime(0  =  0[' 


(17) 


where,   x(t)  =  x(t)  —  x(t)  represents  the  observer's  state 

estimation  error  and  e(t)  =  e(t)  —  e(t)  represents  the  error 

in  uncertainty  estimation. 

We  are  assuming  the  following  for  derivation. 

1.  Matrices  A,  B,  C  are  known.  All  uncertainties  associated 
with  the  matrices  j\  and  B  are  lumped  into  e(t)  . 

2.  The  pair  (A,  B)  is  controllable. 

3.  The  pair  (A,  C)  is  observable. 


KiT. 


OWrver  *■ 


Figure  1 .  Generalized  block  diagram  of  observer  based  MPC  system 

From  (1)  and  (16),  the  dynamic  of  the  observer  error  becomes 
x(t  +  l)  =  (A-LC)x(t)  +  Be(t).  (18) 

which  together  with  (16)  results  in; 


x(t  +  l) 
e(t  +  l) 


A-LC     B 
-MC      0 


x(t) 
e(t) 


(19) 


Simplifying 

(13) 

~x(t  +  l)~ 
e(t  +  l) 

'A    B~ 
0     0 

- 

'  L 
M 

[c  of 

J 

~x(t) 
e(t) 

(20) 


It  follows  that  x(t)  is  an  estimate  of  x(t)  if  and  only  if, 

]j^  X  I R \  <  0  for  q  =  1,  •  •  • ,  (n  +  p)  ■  The  generalized   block 

diagram  of  the  observer  based  MPC  system  is  shown  in 
Fig.  1. 

V.  Mathematical  Modelling  of  Dc  Servomotor 

For  mathematical  model  derivation  of  DC  servomotor,  the 
electrical  &  torque  characteristic  equations  are  formulated. 


—  o 


Figure  2.  Schematic  diagram  of  DC  Servomotor 

Fig.  2,  shows  the  schematic  diagram  of  armature  controlled 
DC  servomotor.  Using  Kirchhoff 's  voltage  law  (KVL), 


V    =RJi+L 


ydtj 


i...  +  k..n._ 


(21) 


where,  km  is  the  back-emf  constant,  R    is  the  motor  winding 

resistance  and  L     is  the  motor  winding  inductance,  /     is 
the  current  through  the  motor  winding. 

Table  I.  Nominal  DC  Servomotor  Parameters  [2 1  ] 


Description 

Symbol 

Value 

Unit 

Motor  electric  time  constant 

*e 

7.74  *  10 5 

5 

Moment  of  inertia 

J, 

2.21*105 

Kg.m 

Motor  maximum  velocity 

u 

max 

298.8 

rad/s 

Motor  maximum  current 

max 

1.42 

A 

Motor  maximum  torque 

T 

max 

0.068 

N.m 

Motor  torque  constant 

K 

0.0502 

N.m/A 

Motor  armature  resistance 

K 

10.6 

Ohm 

Motor  armature  inductance 

Lm 

0.82 

mH 

Open  loop  time  constant 

T 

0.0929 

s 

Open  loop  steady  state  gain 

K 

19.9 

md/(V.s) 

©2013  ACEEE 
DOL01.UCSI.4.1.1064 


26 


vc  ACEEE 


Full  Paper 


ACEEE  Int.  J.  on  Control  System  and  Instrumentation,  Vol.  4,  No.  1,  Feb  2013 


Representing  (21)  in  Laplace  domain  as  follows: 

RJm  +  Lmims  =  va-kana.  (22) 

The  armature  electrical  time  constant  of  system  is  represented 

as. 


R 


(23) 


Expressing  the  torque  equation  into  the  Laplace  domain, 

JeAnS=kJ,n+Tr  (24) 

Combining  (21)  and  (23), 


^Av^  +  ^r5 — _^-^ +  Td(s).       (25) 


eq       m 


R, 


R, 


Transfer  function  from  motor  shaft  position  to  the  motor  input 
voltage  can  be  given  as, 


Gm(s) 


e(s) 


The  plant  system  parameters  can  be  expressed  as, 

Gm(s)=  k" 


(26) 


k  2V 

eq  n 


R 


(27) 


Simplifying  the  above  equation  further, 

r   t\  b  K 

GJs) 


(28) 


where, 


2        and 


a  =■ 


/   R 

eq     m 


J.,,R„ 


k  is  the  steady  state  gain  and  z  is  the  time  constant  of  the 

system. 

Converting  the  transfer  function  into  the  state  space  form, 


0         1 

o   -V 


,B 


0 
K/ 


,c  =  [i  oJd  =  [oJ 


(29) 


VI.  Matlab  Simulation 


A.  QP  Solver 

The  proposed  algorithm  is  coded  in  MATLAB.  Fig.  3, 
shows  per  iteration  computation  time  comparison  of  proposed 
method  with  QR  factorization  and  MATLAB's  QP  solver 
quadprog. 

B.  MPC  Implimentation 

The  designed  MPC  with  PIO  is  tested  through 
experimentation  for  position  control  application  of  a  DC 
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0.06 
0.05 
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1 
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□ h-jjjf-t*.  --Q-^lf  —  a—  <a—  -o — a 0 — 0 —  »  —  ■■■'-■              »—■— ,H 

a — □ — a — «■  -H-  —  o  —  ■** — ■-  *i—  -a — «-■-»  —  »--  h  —  *             '-m-  x^  ^r 

6 


14      1G      n      20 


E       10      12 

Nc.  of  Iterations 
Figure  3.  Per  iteration  computation  time  comparison  of  QP  solving 
methods 

servomotor. 

Substituting  the  parameter  values  from  Table  I  in  (29),  the 
state  space  representation  of  DC  servomotor  model  can  be 
obtained  as, 


A  = 


"0 

i  ■ 

.«= 

~  0 

L° 

-i 

210 

,C  =  [l    0],D  =  [0j 


(30) 


QET  DC  servomotor  model  as  given  in  (30)  was  considered 
for  the  MATLAB  simulation  of  MPC.  The  simulation  results 
of  the  proposed  Observer  based  MPC  and  high  pass  filter 
(HPF)  based  MPC  are  shown  in  Fig.  4,  Fig.  5,  and  Fig.  6. 
1.5 1 p 1 1 < 1 


■"  Reference 
--HP? 


1.5         2         2.5 

Time  (sect 


3.5 


Figure  4.  Comparative  plot  of  position  response 

From  the  results  it  is  observed  that  observer  does  better 
state  estimation  compared  to  HPF  in  transient  time.  Due  to 
better  state  estimation  of  actual  states,  control  input  is 
conservative  and  states  have  less  oscillation  in  the  transient 
response. 
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Figure  5.  Comparative  plot  of  second  state 
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Figure  6.  Comparative  plot  of  control  input 

VII.  Hardware  Implementation 

Efficacy  of  the  proposed  approach  is  demonstrated 
practically  on  Quanser  QET  DC  servomotor  plant  with  0.01s 
sample  time. 
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Figure  7.  Real-time  MPC  Implementation  strategy  for  position 
control  of  QET  DC  servomotor 

The  experimental  setup  has  DC  servomotor  with  Q2  USB 
based  data  acquisition  card  and  a  PC  equipped  with 
proprietary  QuaRc  2. 1  software.  QuaRc  provides  hardware- 
in-loop  simulation  environment  [21].  Motor  shaft  position 
0(t)  is  measured  by  an  encoder.  To  provide  state  estimation 
of  plant,  PIO  is  developed.  The  proposed  real-time  MPC 
implementation  strategy  and  the  experimental  setup  are  shown 
in  Fig.  7  and  Fig.  8  respectively. 

The  performance  of  designed  observer  based  MPC 
scheme  is  compared  with  classical  controller  i.e.  PID 
controller.  Experimentation  has  been  carried  out  with  constant 
load  to  DC  servomotor  as  an  external  disturbance.  Fig.  9, 
shows  the  plot  of  position  response  using  PID  controller 
and  MPC.  It  shows  that  proposed  MPC  scheme  has  better 
transient  response  compare  to  classical  controller. 

Plot  of  control  input  of  MPC  and  PID  controller  is  shown 
in  Fig.  10.  Control  input  of  MPC  is  conservative  as  compare 
to  PID,  which  can  be  clearly  observed  from  the  plot. 
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Figure  8.  Quanser  QET  DC  servomotor  setup 
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Figure  9.  Comparative  plot  of  position  response  using  PID 
controller  and  MPC 
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Figure  10.  Comparative  plot  of  control  input  of  PID  controller  and 

MPC 

Fig.  1 1 ,  shows  the  position  responses  of  the  plant  using 
two  different  prediction  horizon  with  the  control  horizon  kept 

constant  at  N    =3.  From  Fig.  12,  it  is  observed  that 
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Figure  1 1 .  Plant  response  to  reference  position  for  different 
prediction  horizon 
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Figure  12.  Control  input  for  different  prediction  horizon 

N    =10  results  in  a  shorter  rise  time  and  a  shorter  settling 

time  in  the  position  response  then  N  =  40  .  Thus  the 
responses  are  faster  while  decreasing  the  Prediction 
Horizon  N  .  Fig.  12,  shows  that  this  improvement  is  achieved 
at  the  cost  of  larger  control  signal. 

Conclusions 

In  this  paper,  LU  factorization  based  linear  solver  for 
solution  of  IPM  is  presented,  which  is  used  for  real-time 
implementation  of  constrained  linear  model  predictive 
controller  for  position  control  of  DC  servomotor.  The  issue 
of  unavailability  of  internal  state  vector  for  implementation 
of  MPC  is  addressed  by  designing  a  PIO.  Simulation  results 
are  presented  to  demonstrate  the  effectiveness  of  PIO. 
Simulation  result  shows  the  superior  performance  of  observer 
to  estimate  the  state  vector  as  well  uncertainty  compared  to 
HPF  From  experimentation  results  it  is  observed  that,  the 
proposed  scheme  solves  QP  problem  efficiently  within  the 
specified  sample  period  and  performs  well  for  tracking  a 
reference  position  of  DC  servomotor  which  is  relatively  a 
fast  dynamic  system.  Performance  of  the  designed  MPC  has 
been  compared  with  PID  control  via  experimental  results.  It  is 
also  observed  that,  the  controller  is  robust  to  handle  parameter 
uncertainties  and  load  disturbances. 
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